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Abstract 

We derive an expression for the second-order gravitational self-force that acts on a self-gravitating 
compact-object moving in a curved background spacetime. First we develop a new method of 
derivation and apply it to the derivation of the first-order gravitational self-force. Here we find 
that our result conforms with the previously derived expression. Next we generalize our method 
and derive a new expression for the second-order gravitational self-force. This study also has a 
practical motivation: The data analysis for the planned gravitational wave detector LISA requires 
construction of waveforms templates for the expected gravitational waves. Calculation of the two 
leading orders of the gravitational self-force will enable one to construct highly accurate waveform 
templates, which are needed for the data analysis of gravitational-waves that are emitted from 
extreme mass-ratio binaries. 
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I. INTRODUCTION 



When a self gravitating compact object moves in curved spacetime it interacts with its 
own gravitational field. This interaction induces a (gravitational) self-force that the object 
exerts on itself, and thereby affects its own motion. The analysis of this phenomenon is 
much simplified by considering a point particle limit, in which the spatial dimensions of 
the object approach zero. In this limit, the motion of the object is described by a timelike 
worldline, and the dynamics of the object is governed by an equation of motion. 

To derive an expression for the self-force, it is useful to employ perturbations methods. 
In this approximation, one assumes that the object self-gravity can be represented by small 
metric perturbations on a given curved background metric. At the leading order of approx- 
imation, the worldline of the object traces a geodesic in the background spacetime (see e.g. 

At higher orders, the metric perturbations produced by the object, give rise to a self 
force that induces a deviation from the geodesic trajectory. The linear metric perturbations 
that are proportional to the mass of the object fi induce a first-order self-force - proportional 
to 3) 3l- At the next order, the second-order metric perturbations that are proportional 
to /i^, give rise to a second-order self- force - proportional to /i^. 

The main difficulty in deriving an expression for the self-force originates from the fact that 
we consider a point particle limit. In this limit, the metric perturbations become singular 
at the location of the particle, and one is required to introduce a regularization method to 
be able to derive the correct (and finite) expression for self-force. Using different methods 
Mino Sasaki and Tanaka 1^, and independently Quinn and Wald P| have derived the general 
expression for the first-order gravitational self-force in a vacuum background spacetime. To 
be able to go beyond the first-order self-force and calculate the trajectory of an object to a 
higher accuracy, one is required to calculate the interaction of the object with its own second- 
order metric perturbations. These perturbations satisfy perturbations equations with highly 
singular source terms. Recent studies have introduced a regularization method that enables 
o.e to solve these e,.atio. a„d con.t.uet the seeond-o.de. .netne perturbations QQ. In 
this article we follow Refs. |^ |5|] and use the second-order metric perturbations to construct 
an expression for the second-order gravitational self-force. 

The study of self-force has attracted much attention in recent years (see e.g. 0] and ref- 
erences therein) mostly because a calculation of the self-force is needed for future detection 
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of gravitational waves with the planned laser interferometer space antenna (LISA) One 
of the most interesting sources for LISA are extreme mass ratio binaries, in which a compact 
object (e.g. a neutron star or a stellar mass black hole) gradually inspirals towards a su- 
permassive black hole with a mass M (e.g., M//i = 10^). Calculating the self-force that the 
compact object exerts on itself enables one to calculate the inspiral trajectory, and thereby 
prepare gravitational waveform templates for the expected gravitational waves. These tem- 
plates are needed for matched-filtering data-analysis techniques. To determine the binary 
parameters using matched-filtering techniques one is often required to prepare gravitational 
waveform templates with a phase error of less than one-cycle over a year of inspiral j^. 
Calculating waveform templates to this accuracy is a challenging task and requires one to 
take into account the interaction of the compact object with its own second-order metric 
perturbations 0,0]. This provides a practical motivation for the study of the second-order 
gravitational self-force. 

Our analysis is based on perturbation theory. We therefore decompose the full spacetime 
metric g^^^^ as follows 

= 9^.v + 5g^, . (1) 

Here g^^ denotes the background metric satisfying the vacuum Einstein's field equations in 
the absence of the compact object, and 5g^y denotes the metric perturbations produced by 
the compact object. Throughout we shall use g^j^ to raise and lower tensorial indices and 
to evaluate covariant derivatives. Recall that the gravitational self-force originates from the 
interaction of the object with its own metric perturbations. Therefore, the dynamics of the 
object should be described by an equation of motion of the following form 

/^-^^ = f^idt^v, Sg^^, z{t)] , (2) 

where /'^ denotes the desired gravitational self- force ^, z{t) denotes the worldline of the ob- 
ject in the background geometry (at the point particle limit), and r denotes the proper time 
with respect to the background geometry. Recall that the perturbations Sg^j^i, are singular 
at the worldline. This singularity suggests that f^[g^u, ^dfiu, -^(t)] should be interpreted as a 



^ Here it is understood that the dependence on z(t) is other then terms of the form C where C is 
a constant. No observation can distinguish between these terms and the term on the left hand side. 
Therefore, these terms should be absorbed within the term on left hand side in a mass redefinition 
procedure. 
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regularization operator, which takes singular quantities in its input and produces a regular 
self-force as an output. To make this statement meaningful (and useful), we assume that 
depends on the values of Sgfj^iy in the vicinity of the worldline. More precisely, we consider the 
hypersurface of constant time S(r) generated by spacelike geodesies normal to the worldline 
at z{t), and assume that the dependence of /'^(t) on the metric perturbations is in fact a 
dependence on 6g^u[^{T)] in the vicinity of z{t). Throughout this article we shall ignore 
corrections to the self-force which originate from finite size and asphericity of the object. 
These corrections could be analyzed separately using other methods (see e.g. 0,0). For 
sufficiently spherical objects these terms do not affect the first two leading orders which 
concern us here ^. 

We formally expand the metric perturbations in powers of fi which gives 

Sg^,u = /i/i^jz. + + 0{^i^) . (3) 

Throughout this paper we keep the dependence on /i explicit. In the context of perturbations 
theory we can separate the calculation of the self-force into several consecutive steps. In the 
first step, we completely ignore the self- force and consider a geodesic worldline zg{t). Using 
zg{t) as a leading order worldline, we can calculate the first-order metric perturbations h^j,^. 
In this context the first-order self- force is calculated using h^^j and zg{t). The corrections to 
zg{'t) due to the first-order self- force are required as an input for the higher orders terms of 
the self-force, for example they are required for the calculation of the second order self-force. 
Formally expanding in powers of /x we obtain 

f = /^Vfi) [9t^iy, V> ^g(t)] + /iV(2) [9/11^, V' l/^i^^ ^(^)] + . (4) 

Here fl[^ and /^2) denote the first-order self-force and the second-order self-force, respectively. 

To calculate the self-force we shall use a decomposition technique as a part of our deriva- 
tion method (Several authors have previously used a similar technique together with other 
methods of derivation in calculating various types of self-forces, see for example 



m 



mm)- 



Roughly speaking, in this technique one first decomposes the metric perturbations into a 
certain regular piece and a certain singular piece. Next, one calculates the self- force and 



^ In a simpler case of an electromagnetic self-force in flat spacetime it was recently demonstrated that the 
self-force is universal - i.e. at the point particle limit it is independent of the object charge distribution 
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shows (using a certain method) that the self-force is completely determined by the regular 
piece and does not depend on the singular piece. To apply such a technique to the derivation 
of f^^-^ one might attempt to introduce the following decomposition h^^^, = h^^^ + k^p , where 
/i^*"^ is a certain singular piece and K^^^ is a certain regular piece. Here there is a difficulty, 
since this decomposition is nonunique. We can easily generate a family of such decomposi- 
tions by invoking a transformation of the form /i^*"^ h^^^u^ + Q/iu , h'^^u ~^ ^^i^u ~ Qfj-i^ 5 where 
q^u is an arbitrary regular field. Therefore, there is a concern that our final expression for 
fj^^-^ will depend on the specific decomposition in use. In this article we employ a specific 
decomposition (see below) that suffers from a similar nonuniqueness problem. Nevertheless, 
we show that this nonuniqueness does not affect the formulas that we derive for /^'^^ and 

This article is organized as follows: In Sec. |n]we derive an expression for the first-order 
gravitational self-force in a vacuum background spacetime. Tackling this simpler problem 
first allows us to introduce the main principles of our derivation method. In addition, this 
analysis also serves as a check of our method, since we are able to compare our formula 
for /^'^-j with the previously derived expression for the first order self-force P, O]. Next, in 
Sec. mil we generalize the principles which have been introduced in Sec. |nj and derive an 
expression for the second-order gravitational self-force in a vacuum background spacetime. 



II. FIRST ORDER SELF-FORCE 



We consider the first-order metric perturbations in the Lorenz gauge 

where overbar denotes the trace-reversal operator defined by hf^i, = h^^ — {l/2)g^^h^. In 
the Lorenz gauge the first-order metric perturbations satisfy the following wave equation 

□V + 2fly>,, = - 16irr„„ , (5) 

where □ . «-V,V„ denotes the Riemann te,w of the background geo^try 

with the sign convention of reference lL&\. Throughout this article we use the signature 
(— , +, +, +) and geometrized units G = c = 1. In a local neighborhood of the worldline the 
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energy-momentum tensor (at the point particle limit) takes the form of 

/"OO 

Tf,^{x) = fi / gf'{x\zG)uagf{x\zG)u(36^{x - ZG)[-g{zG)r^^^dT . 



OO 



where the four velocity is denoted u'^ = for an arbitrary point x' the notation 



dr ' 

g^" = g^{x\x') denotes the geodetic parallel propagator from x' to x. This bi-vector trans- 
ports an arbitrary vector Aa{x') to a vector = g {x\x') Aa{x') by a parallel propaga- 
tion of the vector along the geodesic connecting x' and x (see Appendix A); g denotes 
the determinant of the background metric, and 6^{x — x') denotes the four- dimensional 
(coordinate) Dirac delta-function. 

Eq. (jH) implies that /^'^^ depends on the metric perturbations h^^, where hereafter h^y 
shall denote the retarded solution of Eq. (0). In a local neighborhood of the worldline 
it is useful expand the perturbations h^^ using the distance from the worldline as a small 
parameter. For definiteness let us consider an arbitrary point on the world line z = ZG{f), 
and expand /i^jy on S(t) as follows 



h.,{x) = -g^g; [c^r^^e-' + + + 0{e')\ . (6) 

Here the overhat indices refer to z, e denotes the length of the geodesic connecting z and 
X. The coefficients jc*'"'- 1 transform as tensors under a coordinate transformation at z, and 
transform as scalars under a coordinate transformation at x. By definition these coefficients 
are independent of e. The explicit expressions for the coefficients {c^.*^^} are not required in 
this paper, but we should point out that some of the coefficients are nonlocal (in particular 
c'^^p i.e., they depend on the entire past history of the worldline (see e.g. jgl). Notice that the 
perturbations h^i, diverge as and therefore the calculation of the self-force f^^-^ requires 
one to introduce a regularization method. 

Our regularization method is based on a decomposition of hf^^ into the following tensor 
fields in a local neighborhood of zg{t): 

V = + C • (7) 

Here /i^j, denotes an instantaneous piece and h^^ denotes a sufficiently regular piece. We 
define an instantaneous piece hj^^ as a piece that admits a local expansion on S('r), reading 



hUx) = -g^gf [4"/^^-^ + <j + 4>] , (8) 



where the expansion coefficients {d^^^\ d^^^, ^i^j} constructed only from combinations of 
quantities in following list: the metric tensor g^xp and tensors constructed from the metric 
and its derivatives (e.g. Riemann tensor and its covariant derivatives), hereafter we shall 
refer to these tensors as background tensors; the four velocity u^; the unit tangent vector 
Va£ which is tangent to the geodesic connecting z and x (see Eq. I A3 1 in Appendix A), 
this quantity transforms as a vector at z and as a scalar at x; and numerical coefficients. 
We define a sufficiently regular piece h^^ as a piece whose ffist-order derivative V ah^J} 
is continuous in a local neighborhood of zg-, and its higher-order derivatives are not too 
singular (More precisely we demand that 6"^"^^ s^h^J}.^, where n > 1, remains bounded 
as X — s> z.). 

Roughly speaking, the instantaneous piece captures the singularity content of h^y, and 
therefore the piece is sufficiently regular. Another important property of hl^p{x) is 
the "instantaneous" property: unlike h^^, the instantaneous piece h^^p^x) on S('r) does not 
depend on the past history of zq, it depends only on instantaneous quantities defined on 
S(f). 

A priori it is not clear that h^^ can be decomposed into an instantaneous piece and a 
sufficiently regular piece. For now let us suppose without proof that this decomposition 
exists. Later (in Sec. HE) we shall provide a specific prescription for its construction. 
Notice that decomposition ((7j) is nonunique, and we are at liberty to introduce an alternative 
decomposition that satisfies the same conditions, reading 

h^.u = {hi, + V) + (^^^ - Ku) . (9) 

where k^^, is both sufficiently regular and instantaneous (i.e. it satisfies all the conditions 
that these fields satisfy), but is otherwise completely arbitrary. At this stage we do not 
provide a specific prescription for the construction of h^^, and h^^}, and our derivation will 
be based only on the general properties that define this decomposition. 

Substituting Eq. (jlj) into Eq. © and using Eq. ((Tj) we find that at the ffist-order 
[ignoring 0{fi^) corrections] the equation of motion reads 

= f^'fii) [9,^^ + ^o{t)] , (10) 

The left hand side of this equation is dimensionless, and therefore the vector /^'^^ must have 
dimensions of {Length)~'^ . Notice that the ffist-order self- force can have a complicated 
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nonlinear dependence on the background metric. However, it must be linear in h^,y, since 
terms which are nonlinear in h^iy must also come with higher powers of fi and therefore belong 
to higher order corrections to f^. Due to this linearity in h^^, we may group the terms in the 
expression for f^^-^ into a piece originating from {Qhu, hj^^u^ ^g{'^)} denoted /(^)/, and a piece 
originating from {gf^^, h^u, ^g(t)} denoted /[J)^^- We formally write this decomposition as 

fii)[9,u, + hl^, zg{t)] = fl{^j[g,u, hi,, zg{t)] + fll^sRia,'^^ ^oir)] . (11) 
Here both /(^^^^ and /^'^^ are well defined vector fields with dimensions of {Length)'"^ . 

A. Instantaneous piece 

We now introduce the principles of our derivation method that is used extensively in this 
article. Let us focus on the instantaneous piece f({-^j[gfj.u,hjj^,,ZG{T)]. For definiteness, and 
without loss of generality, we calculate f^-^^-^j (i.e. fl[^j at z). To derive an expression for /^'^^j 
we first make a list of all the possible tensors (at z) that may be included in the expression 
for We refer to these tensors as the tensorial constituents of /(^)/- In practice, our list of 
tensorial constituents is constructed by studying the possible tensors that can be constructed 
from the quantities {(7^^, h^,, zg{t)}. Next, we combine these tensorial constituents to find 
all the possible vector expressions for /(^)/- In this method we first make a list of individual 
tensors (without their combinations), and consider the combinations only at the next step, 
when we construct the possible vector expressions. Notice that these vector expressions 
must be well defined and must have dimensions of {Length)~'^ . 

The explicit dependence of /^'^^j on zg{t) and g^,, implies that the four velocity and the 
background tensors must be included in the list of the tensorial constituents of fj^-^j- Since 
zg{t) is a geodesic, higher order covariant time derivatives of the four- velocity vanish. The 
dependence of f^-^^j on h^, implies that tensors that can be constructed from hjj^, should also 
be included in our list of tensorial constituents. All these tensors can be expressed using the 
quantities that are listed in Eq. (jHl). We now show that these tensors do not extend our list 
of tensorial constituents. From its definition hj^, depends on both the four-velocity and the 
background tensors. These tensors have already been included in our list. In addition, hj^^, 
also depends on Vo-e. However, the tangent vector Vo-e is well defined only for a specific 
value of X but it becomes multivalued as we change x. Since the self-force at z is independent 
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of x, it can not depend on Va^. Moreover, we can not generate a new well-defined vector 
by taking the limit e ^ of Vq£:. This limit is ill-defined since its value depends on the 
direction from which the limit is taken. We therefore conclude that V aS must be excluded 
from our list. A similar argument shows that e must also be excluded from our list. In 
addition, h^^^ also depends on the parallel propagator g^. Since is a bi- vector rather 
then a tensor, it can not be included in our list of tensorial constituents. However, one can 
construct a well-defined tensor from g^, for example by taking the limit e — > of g^ (or 
by taking the limit e Q oi its covariant derivatives) and thereby construct a tensor at 
z. Notice, however, that g^^ is completely determined from the background metric and the 
geodesic connecting z and x. Therefore, any well defined tensor (independent of x) that is 
constructed from g^ must be a background tensor. We therefore conclude that the tensors 
that are constructed from /i^^ do not extend our list of tensorial constituents. Finally, we 
need to consider the possibility that the Levi-Civita tensor might appear in our list. Notice 
that a term with an even number of Levi-Civita tensors does not introduce any new tensor to 
our list, since these combinations can be expressed using combinations of Kronecker deltas^ 
. It is therefore sufficient to consider a single Levi-Civita tensor in our final expression for 

J(i)r 

Summarizing the above analysis, we have found that the expression for /^'^^^ must be 
composed from the following tensorial constituents: the four velocity, the background ten- 
sors, and a single Levi-Civita tensor. The fact that f^-^^^ has a dimension of [Length)'"^ 
further constrains our list of tensorial constituents. A background tensor with a dimension 
of {Length)'"^ may contain only the following quantities: ^f^j), (y'/ip,^, and f?^^^^- Clearly a 
combination of gfj,o together with gf^y^a can not produce such a tensor, since we may always 
employ a locally flat coordinate system (at z) in which g(j,i,^a = 0. Consequently, a back- 
ground tensor that is constructed from gf^^, g/io^a and has a dimensionality of [Length)'"^ 
must vanish. Therefore, the desired tensor must be linear in q-- - a in addition to a combi- 

provides a uniqueness theorem that ensures us that the 
only background tensor that satisfies all these requirements is the Riemann tensor R^^^^ (up 
to a combination with gfic)- Tensors that include higher-order derivatives of the metric, or 

^ Any two Levi-Civita tensors s^^^^Efj^^prf can be expressed as a sum over terms of the form ±6'^S^6JSfj with 
appropriate index permutations. This sum is often expressed as a 4 x 4 determinant of Kronecker deltas, 
see e.g. jirj. 
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tensors which are nonhnear in the second derivative of the metric [e.g. tensors of the form 
^e-^a/37(5' tensors that are quadratic in the Riemann tensor] have a "wrong" dimensional- 
ity of (Length)'^ ^'^^^^ where n > 0, and are therefore excluded from our list"^ . We conclude 
that each term in fj^^^j must be constructed from a single Riemann tensor R^0a^§ and may 
also include the following dimensionless tensors: u'^, ^f^p, and a single Levi-Civita tensor 

We now show that in vacuum it is impossible to construct a vector from our list of 
tensorial constituents. First we consider vectors that are constructed from: -R^^^^, gjio, and 
u^. We find that all these vectors vanish by virtue of the symmetries Rg^p^^s) ~ 0; ^{a(5)^5 ~ 
and the fact that R^^ = 0. Considering next vector expressions involving the tensor e^^'^'? 
leads to the same vanishing result: An attempt to contract three or four indices of R^^A^g 
with e'^'^P''^ vanishes due to the identity R"'^^^^ = 0, contracting a pair of indices of R^/s^s 
with e^'^^'^ produces a tensor with two pairs of indices (a pair from each tensor). Each of 
these pairs is antisymmetric, and therefore it is impossible to combine the resulting tensor 
with the four velocity and construct a vector. The remaining possibilities vanish trivially. 

Since one can not construct a non- vanishing vector from the above list of tensorial con- 
stituents we are led to the conclusion that 

/fi)/ = . (12) 



Ori have used a similar argument to show that the expression for fj^^^ does not consists of 
any local terms In this article we shall extend this argument in the analysis of f^^) 

(see below). After this analysis was completed we have learned that Anderson, Flanagan 
and Ottewill have further extended and simplified Ori's argument in their calculation of a 
quasi-local expansion of the first-order self-force 



B. Sufficiently regular piece 

Following the method of Sec. IIA we begin our analysis by studying the tensorial con- 
stituents of the sufficiently regular piece of the self-force fj^i^sRldfj.^^ ^^u^ ^ciT)]- Recall that 
fj^i-^gji must be a well-defined vector with a dimension of (Length)^'^ . The dependence of 

^ Notice that tensors with dimensions of {Length)", where n > [e.g. {Rap^sR"'^^^)^'^ Rfii>vp] rnust be 
excluded since they give rise to terms that are ill defined in flat spacetime. 
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f{i)SR "-"^ 91^^^ ^^'^ ^g{'^) implies that the background tensors and the four-velocity must be 
included in our list of the tensorial constituents of f([-^sR- However, in Sec. IIA we found 
that in vacuum these tensorial constituents (together with the Levi-Civita tensor) can not 
be combined to give a non- vanishing vector with a dimension of (Length)'"^ . Therefore, all 
the terms in the expression for /(^^^^ must contain an explicit (linear) dependence on h^^. 

Since h^l^ is sufficiently regular it admits the following local expansion on S('r) (see 
Appendix A) 

h'Jix) = -gj^-gf [hll - hlf.ee-'-' + 0{e')] . (13) 

For illustration consider the simplest case of static particle in flat spacetime. In this case 
Eq. (|13|) is reduced to a Taylor expansion: using Cartesian coordinates = {x^,x^,x^) we 
flnd that in this case we have g^ = 6J^, e = (a;" — z")(x^ — z'')Sab, and ee''^ = —{x'^ — z'^). 

We now examine the terms in Eq. (fT^ that may appear in our list of the tensorial 
constituents of ffi-^g^- Following Sec. IIA we flnd that the quantities g^^ , e, and e''^ that 
appear in Eq. ()13|1 are excluded from our list. We are therefore left with the terms h'^^, 
^ff-y' possibly higher-order derivatives (if they exist). Recall that in our notation 
h'^^ has a dimension of (Length)'^ [see Eq. (jS))]. Therefore, higher-order derivatives of 
the form h^f^g g , n > 1, must have dimensions of {Length)~'^~"' . By dimensionality, 
these terms must be excluded from our list 21]. Consequently, from the various terms 
that are listed in Eq. |T3|l only the terms h^^ and h^^^^ can appear in the expression 
for f([^gji- For example we can construct the following vectors {g""" + u°'u'')u'^u^h'p^?.^ and 
{g"^ + u°'u^)u'^u^hp^?.~. Since the set of all possible vector expressions that may contribute to 
f{i)SR empty set, we can not immediately implement the method of Sec. IIA. To 

get around this difficulty we employ a gauge transformation. Previous studies have shown 



221 and also discussion in section 



that the self-force is a gauge dependent quantity (see 
lie below). Therefore, we may employ a gauge transformation to simplify our analysis. 
Consider a regular gauge transformation of the form x^ ^ x^ — We denote the first- 
order perturbations in the new gauge with /i^^. In this gauge we have h^^ = h^^^ + h^!}^^\ 
where h^u^^^ = h^Jj- + + ^j,;^. Here we have included the entire gauge transformation in 
the definition of the new sufficiently regular piece h^u^^\ thus leaving hj^^^ gauge invariant. 
We demand that C,^;up "will be continuous in a local neighborhood of the worldline, and 
thereby we guarantee that h^!}^^^ will satisfy the conditions of a sufficiently regular piece. 
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We shall refer to the new gauge as Fermi-gauge. In Fermi gauge Eq. reads 

nii)\9,u. hi + hl^'^^\ ZG{r)] = fl^,[g,^, hi^, zo{t)] + fl^sRiF)[9,u, C^^^, zcir)] , (14) 

where fp^^^^ is the first-order self-force in Fermi gauge and fl{)sR(F) corresponding 
sufficiently regular piece. We impose the following gauge conditions along the world line 

[C^^^] - + e.;. + J = , (15) 

- [K^ + e.;^ + ^^;J;J = • (16) 

In complete analogy with the Lorenz gauge, each term in f(\^sji(^p) (in Fermi gauge) must 
explicitly depend on one of the first two coefficients in the expansion of h^u^^^ [similar to 
Eq. (fTSj) ]. Since these coefficients vanish all these terms must be equal to zero. We therefore 
conclude that 

f{l)SR(F) = • 

We now briefly explain how to construct Fermi-gauge. By contracting Eq. (|16|) with 
u^, we flnd that -^h^u^^^ = 0. This result is consistent with Eq. (|TH|l . We can now solve 
Eq. (fT3j) at an initial point zg{tq). Satisfying Eq. (fT^ guarantees the validity of Eq. (fT^ 
everywhere along Zg{t). We construct an arbitrary gauge vector ^(o)^ at Zg{tq), and demand 
that 

C(o)fi;u = ^iQ)u■,^l = --h'^ui'To) , (18) 

thereby satisfying Eq. (fTH|l . We now turn to Eq. (fTT)|) . Viewing these equations as algebraic 
equations, we flnd that we have a set of 40 equations for the 64 variables C,fj.;up- By introducing 
the commutation relation 

2^^;[H = R\uM ) (19) 

we include another 24 equations which brings us to 64 equations, as desired. Using the 
identities of the Riemann tensor we obtain from Eqs. p6|19j) the following relation 

^v\iia = R afiu^^ ~ ^^i/fia ' (20) 



^ Fermi gauge was first introduced in Ref. In this reference the gauge vector was defined using Eqs. 
P5llt)|) with the replacement h^^ /i^^, where hjj]^ is Detweiler and Whiting regular piece 0. These 
different definitions coincide once we express h^^} in terms of h^[^ in Eq. below. 
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where 

^ySR = Ifh^R +h^R -h^R \ 

Here all quantities are evaluated on the worldline. By contracting Eq. (j20|) with u°'u'^ we 
obtain 

L = R^.Mu'-u^ - w-w'^^rf^ , (21) 

where = -^y^- One can solve this second-order transport equation by using ^(o)/^, and ^(o);^ 
as initial conditions, and thereby obtain a gauge vector .^^(t) along the worldline. Similarly 
we can construct a first-order transport equation for ^^.j^ by contracting Eq. ()20|) with m". 
Using i^ij) together with Eq. (|TH|l this first-order transport equation can be integrated to 
give i^-^ui,''')- iii.Ma{T) is then obtained by substituting ^/^(r) into Eq. ((201) • This completes 
the formal solution of Eqs. ()15pi(i|l . and allows one to calculate the leading terms in a local 
expansion of in the vicinity of the worldline. This local expansion could be arbitrarily 
continued, thereby obtaining a global definition for Fermi gauge. 

C. First-order self- force in Fermi gauge 

Combining Eqs. ()12ll7p together with Eq. p4p we find that in Fermi gauge the first-order 
self-force vanishes, namely 

/^(i) = . (22) 

This result might appear surprising at first. The surprise is removed once one recalls that 
within the framework of general relativity the complete information about the motion of a 
point like object includes two pieces of information: the object's worldline together with the 
spacetime metric. It is therefore not surprising that one may transfer information about the 
motion of the object from the description of the worldline to the description of the gauge. 
This also implies that an expression for the self-force is physically meaningful only if the 
gauge is specified. 

D. First-order self-force in the Lorenz gauge 

We now derive the expression for the first-order self- force in the Lorenz gauge. For this 
purpose, we consider the inverse gauge transformation i.e. a transformation from Fermi 
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gauge to the original Lorenz gauge, generated by a;" + /i^", where satisfies Eqs. 

()15ll6p . Barack and Ori have derived the gauge transformation formula for the first-order 



self-force 



22l |. Their analysis shows that under a regular gauge transformation of the form 



— >■ x° — /i^" the first-order self-force transforms as 

fii) ^ /« + ^fm ' (23) 

where 

Sfii) = -(^?"' + «"^')eA - R'^.x^.u^en' . (24) 
Using Eqs. ()22I23I24|) we find that in the Lorenz gauge the self-force is given by 

/« = fm - 5fii) = -5fw ' (25) 

where the minus sign originate from the fact that we are considering the inverse gauge 
transformation. We now calculate /^"^ by first substituting Eq. (j2H) into Eq. (|^. and 
then substituting into Eq. (j2Sl)- In this way we obtain the following expression for the 
first-order self-force in the Lorenz gauge 

/i^S) = f^'K^'-'^h^J;, . (26) 

where 

We have already mentioned that decomposition ((7j) is nonunique, and one can generate a 
family of alternative decompositions 0, where k^y is both sufficiently regular and instanta- 
neous, but is otherwise completely arbitrary. Clearly our result in Eq. fl26p is self-consistent 
only if it is invariant under the transformation —>■ h^^—kap. This invariance follows from 
properties of kap. Since fco/j is sufficiently regular the vector K^°'^'^ k^p.^ is well defined on the 
worldline, and since kap is instantaneous it follows from subsection IIA that K^°'l^^kai3-^^ = 0, 
meaning that Eq. is invariant under the transformation h^j^ — > h^p — /c^^. 

E. Specific construction of the metric decomposition 

So far we did not specify how to practically construct decomposition ((7j), nor did we show 
that this decomposition exists. In this subsection we complete our derivation by providing 



14 



a specific prescription for the construction of decomposition ((Tj). Decomposition ((Tj) can be 
constructed in more than one method. For simphcity we choose to follow a previous analysis 
by Detweiler and Whiting [3| and relate decomposition ((Tj) to their decomposition, reading 



hap = h% + h^p. (27) 

Here the (trace reversed) singular piece /if^ satisfies the inhomogeneous wave-equation 
In a local neighborhood of the worldline it is defined by 

/OD 
G'^^^p[x\ziT)]u'^{T)u^{T)dT, (28) 
-oo 

where the singular Green's function is given by 

G^ua'A^l^'] = l[u^ua'p'[x\x']6{a) + Vf,^o.> p>[x\x']e{-a) . (29) 

Here a = a{x\x') denotes half the square of the invariant distance measured along a geodesic 
connecting x and x'; U^^a'f3'[x\x'] and V^^a'i3'[x\x'] are certain regular bi-tensors (for their 
definitions and properties see e.g. 0, [2^), denotes a step function. Eq. implies that 
for a fixed x only a finite interval contributes to the integral in Eq. (f^ . This interval 
approaches zero as as the evaluation point x approaches point z on the worldline. The 
regular field /if^ satisfies a homogeneous wave equation [Eq. (0) with T^^, = 0], and is a 
smooth field in the vicinity of the worldline. 

Expanding h^^^ on the hypersurface of constant time S('r) gives (see Appendix B) 

Comparing Eq. (jHIH) with Eq. (jHl) we find that the leading terms in the expansion h^p{x) 
[and h^p{x)] satisfy the conditions of an instantaneous piece. We therefore identify h^^pi^) 
with the expression in Eq. (jHUj) up to 0{e) inclusive. This identification implies that in the 
vicinity of the worldline we have 

hi^,=h% + 0{e'), (31) 

where the 0{e'^) discrepancy between and /i^^ originate from terms of order which 
appear in h^^ but are absent from h^p. Eq. (j3T|) implies that the expression that have 
been identified as satisfies the required condition of a sufficiently regular piece^. We 



^ Notice that the singularities of the derivatives of the 0{£^) term in Eq. H31() are constrained by the form 
of the O(e^) terms in Eq. ljB4p in Appendix B, as required. 
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therefore conclude that we have constructed a decomposition which satisfies the conditions 
of decomposition ((7j). Substituting Eq. into Eq. yields the standard expression for 
the first-order self-force 



(32) 



F. Summary 

Before we generalize our method and tackle the second-order problem it is instructive to 
summarize the derivation method that we have used. First, in Eq. (0) we introduced the 
general form of the equation of motion. Next, in Eq. ((7j) we introduced a decomposition of 
the first-order metric perturbations into instantaneous and sufficiently regular pieces, and 
defined their properties. Using these properties we showed in Sec. IIA that one can not 
construct a well defined vector (with the desired dimensionality and the correct scaling in 
/i) from the tensorial constituents of the instantaneous piece. Next, in Sec. IIB, we invoked 
a gauge transformation from Lorenz-gauge to Fermi-gauge, and used the same method (as 
in Sec. IIA) to show that the sufficiently regular piece of the first-order self-force vanishes. 
At this stage we concluded that in Fermi-gauge the first-order self- force must vanish. Next, 
in Sec. IID, we used the inverse gauge transformation together with Barack and Ori trans- 
formation formula to derive the expression for the first-order self-force in the original Lorenz 
gauge. Finally, in Sec. HE, we employed Detweiler and Whiting decomposition and provided 
a prescription for the construction of decomposition ((Tj). Using this prescription we have 
obtained the standard expression for the first-order self-force in the Lorenz gauge, given by 
Eq. ()32j) above. The fact that the entire derivation followed from the general properties 
of decomposition ((Tj) led to a simple analysis. This simplicity enables us now to tackle the 
more complicated second-order problem. 



III. SECOND ORDER SELF-FORCE 



From its definition the second-order self-force depends on the second-order metric per- 
turbations. It is therefore necessary to begin our discussion by introducing a construction 
method for these second-order perturbations. Already this preliminary stage introduces dif- 
ficulties. The difficulties and their resolution were studied in Refs. 0, |^. We now briefiy 
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summarize the main results of these references. 



A. Construction of the second-order metric perturbations 

Following Refs. [J, |^ we specialize to a compact object which is a Schwarzschild black- 
hole with a mass /i. We focus on the region far from the black hole i.e. at distances larger 
than ve, where ^ fi, and refer to this region as the external zone. In the external zone 

the geometry is dominated by the background geometry, and the full metric is represented 
by Eq. ([T]), where the metric perturbations are given by Eq. Q. We let r£;(yu) — as 
/i ^ 0. In this limit we write the perturbations equations as 

D^,[h]=0 , x^^, (33) 
D,,\[] = S,,[h] , x^j. (34) 

Here 7 is a timelike worldline, the operators Df^,y and 5*^,^ are obtained from an expansion of 
the full Ricci tensor, where for brevity we have omitted tensorial indices inside the squared 
brackets. This expansion is obtained by substituting Eq. (P) into the Ricci tensor, which 
gives -R^"" = -R^fJ + -R|j1] [Sg] + R^w [Sg] + 0{6g^) (For explicit expressions of the terms in this 
expansion see e.g. jl^-)- To simplify the notation we denote R^f^}[h] with Df^iy[h], and denote 
—R^^u[h] with S^i,[h], where h^^jy is expressed in terms of hf^iy. D'Eath have shown that in 
the limit, the perturbations h^i, are identical to the retarded first-order metric-perturbations 
produced by a unit-mass point particle tracing the worldline 7 At the leading order of 
approximation, 7 is a geodesic of the background spacetime Q , which we continue to denote 
with zg{t). Therefore, in the Lorenz gauge the perturbations /i^i, are given by the retarded 
solution of Eq. (jSJ. At the next order of approximation, the worldline 7 has an 0(yu) 
acceleration, induced by the first order self-force. Hereafter we shall adopt Fermi gauge for 
the first-order perturbations. Recall that in this gauge the first-order self-force vanishes and 
the geodesic worldline zg{t) becomes sufficiently accurate up to acceleration of order 0(/i), 
inclusive. In this gauge we may replace 7 with Zg{t) when we construct the second-order 
metric perturbations \4,M. 

nn 

Refs. |5|| show that the second-order metric perturbations can be decomposed as follows 

Ifiiy = ^i,u + Slf,^. (35) 
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where the first piece ifj^iy is given by the following linear combination of terms that are 
quadratic in h^j^ 

i'f.u = ^ [2{cA^i + CB<^?J - 7(cc<^J, + cn^f!,)] . (36) 

Here = h^" ,h%, V% = h^%hl, , y,^, = {h^^Ph\^)g,,, y,^, = {h^%f g,.', and the 
constants ca, cb, cc, cd must satisfy ca + cb = 1, cq + cr, = 1, but are otherwise arbitrary. 
The second piece 5Tfj_u is given by the retarded solution of the following wave equation 

U6l^, + 2W;,6^p = -2SS^, . (37) 

Here 5S^,^ = S^y - D^,y[^]. 



B. General considerations 

The second term in Eq. namely 

^^^f{2) [9t^v, ^i^y^ zg{t)] (38) 

provides a complete list of the quantities that the second order self-force may depend upon. 
Notice that working in Fermi gauge justifies the substitution z{t) —>■ zg{t) which was made 
in Eq. (p?H|) . This equation also states that the second-order self-force is proportional to 
/i^, meaning that this self-force induces an acceleration which scales as /i^. To obtain the 
desired fi^ scaling of the acceleration, the second-order self-force must be linear in fi'^lfj,^- 
Higher powers of fi'^lfj,^ are therefore excluded from /^'^^. Similarly, terms in /^'^^ that depend 
explicitly on /i^^ must be quadratic in hj^^. Combining all the quantities that /^2) depends 
upon together with Eq. we obtain 

f{2) = f{2)Ai9ap, ZG{r), h^g] + /g)B[c/a/3, Zcir), 6l^s] ■ (39) 

Here the first term is quadratic in /i^^ and the second term is linear in 6l^s- 

From Eqs. (|2l4p we find that /^'^j has a dimension of {Length)~^ . Following Sec. II, we 
shall list all the tensorial constituents of f^^y The dimensionality of f^^^ restricts our list of 
tensorial constituents to tensors with dimensionality of (Length)'"' , < n < 3. This follows 
from the fact that we do not have tensors with dimensionality of [Lengthy , j > at our 
disposal. 
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C. First term fl^2)A 

First we focus on the term fl'^^^[gai3,ZG{T),h^^] in Eq. Here it proves simpler to 

consider a general dependence (not necessarily quadratic) of /^2)yi '^^ ^f^u- the following 
analysis we shall often use decomposition (j77|) . We write this decomposition in Fermi gauge, 
and obtain 

haf3 = hafS + h^js"^ , h^jf^ = h^p + + , (40) 

where satisfies Eqs. ()15I16|) . Using this decomposition we schematically write the first 
piece as 

We shall now list all the tensorial constituents of fl^2)A- First, we shall consider the quantities 
{9ai3, zg{t), h^^}, list the relevant tensors, and discuss the possible expressions for /(2)a 
that can be constructed from these quantities . Next, we shall consider all the remaining 
possibilities. 



1. terms 

We now list all the tensorial constituents that can be constructed from the quantities 
{9ai3, ZG{j),h^g} and have dimensions of (Length)''^, < n < 3. From the quantities 
{9af3,ZG{T)} we obtain the following hst: u^', gf,o, R&p^s^ and -R^^^^.,- In addition, our 
list includes the Levi-Civita tensor. Recall that /i^^ admits the local expansion The 
coefficients of this expansion up to 0{e) (inclusive) satisfy the conditions of an instantaneous 
piece. More generally, the higher order coefficients in this expansion are composed from the 
same kind of instantaneous quantities, namely: background tensors, m'^, Vq,£:, numerical 
coefficients (see Appendix B). Following the arguments of Sec. IIA we now find that h^^ 
does not add any new tensor to the above list of tensorial constituents. Examining this list 
we find that a vector with the desired dimensionality of {Length)"^ must be linear in R^^^x^^.^ 
and may also include combinations of the following dimensionless tensors: m^, ^f^j) and e^l^^^ , 

Examining our list of tensorial constituents shows that they can not be combined to 
give a vector with a dimension of (Length)^^: Consider first vectors that are constructed 
from the following tensors R^f^^s-^, u'^, and gfio- By virtue of the symmetries R£,0(^x^s) ~ 0' 
R{&0)^s ~ ^ ^^'^^ ^^^^ R&fi ~ 0' fi^*^ t'^at a siS'^ "^^^y possibility is 
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i?*^^^- jlI^u^ . This possibility vanishes by virtue of the Bianchi identities which imply that in 
vacuum we have -R^^^^ - = 0. Similarly, using the properties of the Riemann tensor together 
with properties of the Levi-Civita tensor one finds that all the vector expressions composed 
from the tensor e"^^^ together with i?^^^^ gj^o, and m'^, vanish as well. We conclude that 
by combining tensorial constituents that originate from the list {QaPiU^ ih^fs) "^^^ 
construct a vector with dimensionality of (Length)^^ . 



2. Remaining terms 

Since the tensors {ga/s, u'^, h^s} can not produce a vector with the desired dimensionality 
we find that each term in fl^2)A i^ust explicitly depend on h^jf\ Expanding /if^^"* in a local 
neighborhood of the worldline gives 



2 a/3;7'5 



(42) 



Notice that the leading terms that scale as and vanished by virtue of Fermi gauge. In 
this expansion only the coefficient of the term has a dimension which fits our dimension- 
ality condition of {Length)~" , < n < 3. Therefore, only this coefficient can contribute 
to our list of tensorial constituents. To construct a vector with a dimension of (Length)'^ 
we must combine h ~ - with the other dimensionless tensors from our list, namely m^, Oup, 
and £"'f^^^. One can easily construct such a vector, for example {g""^ + u'^u^)u^u^u'^h^^^g. 
Following Sec. IIB we use the gauge freedom to eliminate all the possibilities of constructing 
such a vector. We extend the definition of Fermi gauge by restricting the values of i^-vapiT) . 
From Eq. (j4Uj) we obtain 

= ^^i/;q/3 + ilJ.\uaP + iu;iial3 ■ (43) 

Following Eq.ipUj) we introduce the anzats 

Here all the quantities are evaluated on the worldline, ST^^^ = ^{h^^y.^ + h^a^^ — /i^^.j,), and 
Du^a/s is a tensor field which is defined on the worldline (see Appendix C) . Substituting Eq. 
into Eq.pHj) gives 

Cfi = 2^(M-)"/3 ■ (45) 
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Unfortunately, D^^oip is not completely arbitrary and its values are constrained by various 
identities [e.g. identity (fT^ ]. These constrains prevent us from annihilating h^l^^p by 
invoking a gauge transformation. In Appendix C, we choose Dy^^p such that all the required 
constraints are satisfied. We shall refer to this specific choice of gauge as extended Fermi 
gauge. In this gauge we have 

where P^J^^ ^ is a certain combination of and 5^ (see Appendix C); and nSR^j^^j is the 
linear term in the following expansion in powers of /x 

Here the Riemann tensors are evaluated using the metric inside the squared brackets. Mean- 
ing that fiSR^^^^ is the linear perturbation of the Riemann tensor, evaluated with the metric 
9fiu + l^h^u^^ Notice that SR^^^^j has the same algebraic symmetries as the Riemann tensor. 

Furthermore, in Appendix C we show that in vacuum is traceless. We now follow 

0/370 

the arguments in Sec. IIA. There we showed that in vacuum no vector can be constructed 
from the tensors -R^^^^, u'^, e"^^^, and Similarly, the above mentioned properties of 
6R^^;^l imply that in vacuum no vector can be constructed from 6Rf[^} , m", £'='i^'y\ and 0-3. 
Employing Eq. ()46|) we now find that in the extended Fermi gauge we have 



D. Second term 1^^2)3 

Here we focus on the second term f\^^-^^[ga(3,ZG{T),5l^s\ in Eq. (jH^ . Following Sec. II, 
we decompose 5l^s^ in a local neighborhood of zg{t), into a generalized instantaneous piece 
and a sufficiently regular piece 51^^, reading 

51,, = 5ll + 5ll^ . (49) 

Here 5/^^, is defined as a piece that admits the following a local expansion on S('r) 

= 9a%' [e% + + ^ logs] , (50) 
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where the expansion coefficients {e^, e\^l, e^j} are composed only from the quantities in the 
following list: background tensors (at z), m^, Vq,£:, ^R^^-^j and numerical coefficients. We 
define a sufficiently regular piece 51^^ such that its first-order covariant derivative V a^l^J} is 
continuous in a local neighborhood of zq and its higher-order derivatives are not too singular 
(More precisely, we demand that e~^^^V 5^51^^.,^, where n > 1, remains bounded as 
X z.). 

We initially assume that decomposition ()49|) exists. Later, in Sec. IIIF, we shall provide a 
specific prescription for its construction. As in Sec. IID, the properties of decomposition ()49|) 
imply that this decomposition in nonunique. Following Eq. (jHI), we can generate a family 
of decompositions by invoking the transformation 5/^^, 61 jj^^ + 61^^ 61^^ — j^u, 
where j^u satisfies the definitions of both 5/^^, and 61^^. Following the discussion in Sec. IID 
reveals that the non-uniqueness of decomposition ()49|) does not affect our final expression 
for the second-order self-force [see Eq. (j^Hjl below]. 

Using Eqs. ()39I48I49|) together with the fact that /^2) niust be linear in the second-order 
perturbations we find that we can decompose /^2) follows 

/(2) = /(2)b[^"/3' ^Cir), 5l^s] = f^2)l[9af3, Zcir), 5ll^s] + f^2)SRi9af3, zg{t), 51^^] . (51) 



1. Instantaneous piece 

First we consider the instantaneous piece f^2)i[9ai3, -2g(t), ^li/sl ^'i- (E))- The quantities 
da/B, zg{t), 5ll^^ in this expression together with the definition of 5/^^ reveal that the tensorial 
constituents of /^2)/ a^^: the background tensors, the four velocity, the Levi-Civita tensor, 
and the tensor 5R~ . Recall that in Sec. IIIC we showed that in vacuum these tensors 
can not be combined to give a vector with the desired dimensionality of {Length)~^ . We 
therefore conclude that 

ft2)i = • (52) 



2. Sufficiently regular piece 

Next we consider the sufficiently regular piece f(^2)SR- Following the analysis in Sec. IIB, 
we list the tensorial constituents of f(2)SR\-9<^P^ ^g(t), ^l^s^]- The explicit dependence of f'(2)SR 
on the quantities {gap., zg{t)} implies that the background tensors and the four velocity must 
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be included in the our list of the tensorial constituents of fl^2)SR- ■'■^ addition, this list must 
also include the Levi-Civita tensor. Recall that in Sec. IIIC we showed that in vacuum these 
tensors can not be combined to give a non- vanishing vector with a dimension of {Length)~^ . 
Therefore, each term in the expression for fl^2)SR ^^^^ depend (linearly) on 61^^. 
Since 61^^ is sufficiently regular it admits the following expansion on S(f) 

= 9a%' [Sllf - Sllf-ee-''^ + 0{e')\ . (53) 

Recall that Sl^^ has a dimension of ( Length)'^. Thus Sl^^^~, ^ n > 1 must have dimensions 
of (Length)"^'"' . By dimensionality, these higher-order derivatives must be excluded from 
our list [21]. Consequently, each term in the expression of f'(^2)SR ^^st depend on 51^^^ 
or To eliminate all the possibilities of constructing a vector with a dimension of 

{Length)~^ we employ a purely second-order gauge transformation of the form x'^ — > — 
/i^^|^2)- The second-order perturbations in this second-order Fermi gauge are given by Ij^^, = 
^ij.u + ^{2)fj.]u + ^(2)i/;^c Similar to Sec. IIB, we include the entire gauge transformation in the 
definition of the new sufficiently regular piece 6l^u^^^ = 61^^ + ^(2)^*;;^ + i{2)v;n and demand 
that ^{2)fi;up will be continuous in a local neighborhood of the worldline and thereby ensure 
that 51^!}^^^ will satisfy the conditions of a sufficiently regular piece. In the second-order 
Fermi gauge Eq.(jHI]) reads 

/f(2) = /(2)/[^"/3' ^g(t), 5ll^s] + f^2)SR{F) [^"/3' ^g{t), 51^^^^'^] . (54) 

where fp^^) is the second-order self-force in the second-order Fermi gauge, and /(2)5r(f) 
corresponding sufficiently regular piece. Following Sec. IIB we impose the following gauge 
conditions along the world line Zg{t) 



[«?'^.»,„=O.K.T'].„„ = 0- (55) 

i,r,sar/ must depend on 61^^ or 6i:'\. 
end on 61^^^^'' or 61^ 

ap a 

virtue of Eq. (fS^ . We therefore conclude that 

f(2)SR{F) = • (^^) 



Recall that each term in fl^2)SR ^^ust depend on 61^?^ or ^l^J^.y In complete analogy, each 
term in ft^aatr^', must depend on 61^^^^^'' or 6l^^^^^\ Notice that these coefficients vanish by 
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E. Second-order self-force 



Combining Eqs. ()52I56|) together with Eq. ()54|) we find that in the second-order Fermi 
gauge the second-order self-force vanishes, namely 

= . (57) 

We now invoke an inverse gauge transformation of the form x'^ — ^ -|- /i^^[^2) calculate 
the second-order self-force in our original second-order gauge. For this purpose we need to 
generalize Barack and Ori gauge transformation formula of Ref. 2^ [see Eqs. fl23l24|) ] to 



second order. The generalization turns out to be trivial. The first-order expression ()24|) is 
obtained by calculating the leading order acceleration of the worldline which is induced by 
a first-order gauge transformation. Here we consider a purely second-order gauge transfor- 
mation of the form — /i^^^^), but similar to Ref. [2^, we are still calculating the 
leading order acceleration which this transformation induces. Therefore, the final expression 
has the same form as Eqs. ()23I24|) . with the following substitutions C,'^^)^ f(i) ~^ /(2)' 

The final piece of the calculation directly follows from the analysis in Sec. IID. We use 
the gauge conditions to derive an expression for ^^'2)) use this expression to find 
Sf(^2) /(2)- This calculation produces our main result which is the following expression 
for the second order self-force, reading 

/.VS) = ■ (58) 

Here the first-order metric perturbations satisfy the conditions of the extended Fermi gauge, 
the second-order metric perturbations are constructed using the prescription of Sec. IIIA, 
and the piece 61 ^^i, is decomposed according to Eq. . To complete the derivation we now 
provide a prescription for the construction of this decomposition. 



F. Specific decomposition of 51^^^ 

To construct decomposition (j49|) we need to study the singular behavior of 5/^,^ in the 
vicinity of the worldline. For this purpose, we first study the singular behavior of the 
source term —26Saf3 of Eq. (jSTj). From its definition 6Sai3 has the schematic form 6S = 
V h^V Sz h^W , where schematically denotes the first-order metric-perturbations 
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in the extended Fermi-gauge, and denotes "and terms of the form Recall that 
diverg GS clS S clS £ 0. Therefore, the schematic form of 6Sai3 suggests that 6Sai3 should 
diverge as However, Refs. 0, 0] show that the leading singular terms in SSap cancel 
out and 6Sa/3 diverges only as e~^. Formally expanding —26Sa/3 on S('r) gives 

-26S^p = gJ^gf[e-'A^^ + e-'B.^ + 0{e')] , (59) 

where the coefficients A^^, B^^ are independent of e. Notice that this equation uniquely 
defines A^^ and B^^. We now list the tensorial constituents of the coefficients A^^, To 
construct this list, we substitute decomposition pUj) into the schematic expression for 6Sa/3 
and use Eq. (jl^ . Eq. ()30|) . and Eq. ()B4|) in Appendix B. We find that each term in the 
expression for A^^ must be linear in Rxfjps-y ^^^^ each term in the expression for B^^ must 
be either linear in SR- or linear in VnRs.---. In addition to these tensors, the coefficients 
A^^ and B^^ include only combinations of the following dimensionless quantities: u^^, gf^y, 
Va£, and numerical coefficients. 

To construct 5/^^ we first construct a solution to the following equation, defined in a local 
neighborhood of the worldline 

U^^^ + 2R\%^,, = -gf-g/ie-'A-^ + e'^B^^ + 0,(e°)] . (60) 

Here Oa(£°) denotes an arbitrary quantity of order e^, meaning that we do not restrict the 
values of the 0{e^) terms on the right hand side. To construct a solution in the desired form 
we substitute 

^al3 = 9a9p^ap (61) 

into Eq. (jUUI) . and solve directly for (^^^. Notice that A^^, 5^^, and (^^^ transform as scalars 
under a coordinate transformation at It is useful to decompose these quantities into scalar 
spherical harmonics 

oo I 

^<i/3 = E E yUe,^)4>%{e.t), (62) 

i=0 ■m=~l 
oo I 

;=o m.=-l 
oo I 

1=0 m=-l 
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To define the variables {t,6,(f), consider Fermi- normal coordinates based on the worldline, 
where (a = 1, 2, 3) denote the spatial coordinates and t denotes the time coordinate. 
The angular variables are defined through the canonical angular parametrization, namely 
= £:(sin^ cos sin^^ sin 99, cos^^). In Appendix D we construct the following solution to 
Eq. (jnOI), reading 

m = Ifooo , / = (63) 

(^i'? = --al"^ + -bl'^eloge ,1 = 1 

05 = -[/(/ + l)]-ia5 + [2 - /(/ + , / > 1 . 

Notice that the angular dependence of A^^ and B^^ is completely described by combinations 
of the quantity VqE. Therefore, when we calculate the spherical-harmonics coefficients 
and 6^, we end up integrating only over combinations of Vq,£. Therefore, (^^^ is composed 
form the same well defined tensorial constituents that appear in the expressions of and 
B^0. In addition, ip^^ also has a non-trivial angular dependence which can be expressed 
using combinations of Va£- Comparing Eqs. ()6ip62|63|l with Eq. (jKHjl we find that i^ap 
satisfies the conditions of a generalized instantaneous piece. We therefore identify (fap with 
the desired generalized instantaneous piece, namely 

Slip = ^ap . (64) 
By subtracting Eq. (plj) (for the above solution i^ap) from Eq. (j?fj) we obtain 

ni6Lp - if^p) + 2R\%{51,, - ifrjp) = ^/^/[0(£°)] . (65) 

In Appendix D we use the fact that the source term of Eq. ()65|) is non-divergent to show 
that 6lap — ffap satisfies the conditions of a sufficiently regular piece. We therefore make the 
following identification 

Sl-aP = SlctP — ^aP , 

which completes the construction of decomposition (j49j) . 
G. Summary 

Eq. ()58p provides a formula for calculating the second-order gravitational self-force in a 
vacuum background spacetime, given the sufficiently regular piece Sl^^. To use this formula 
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one has to be able to calculate both first-order and second-order metric perturbations, in 
specified first and second order gauges. We now briefiy summarize a prescription for the 
construction of Sl^p. 

The first step in this construction is to calculate the retarded solution of Eq. (0) 
which provides us with h^^ - the (traced reversed) first-order metric-perturbations in the 
Lorenz-gauge. The next step is to calculate the gauge vector that generates the gauge- 
transformation h^^ — > hfj_^ + ^y-fj, + from the Lorenz-gauge to the extended Fermi-gauge. 
To be able to calculate one first has to calculate the regular piece h^j^^zc) together with 
its first and second order derivatives at the worldline. This preliminary calculation follows 
from Eqs. (j27|28l29p . Next, one follows the prescription in Sec. IIB [immediately after Eq. 
(fTTj) ] and calculates ^^{zq) together with its first and second order derivatives. To construct 
the third order derivative of C,^ one should use Eq. ()44|) together with Eqs. ()C6IC7|) in 
Appendix C. Once ^^(-Zg) together with its first, second and third order derivatives on the 
worldline are obtained, one can use these quantities to construct a local expansion for ^'^ 
(this expansion can be continued in an arbitrary manner away from the worldline). 

The next step in our construction is to calculate the (traced reversed) second-order metric 
perturbations Ifj,^. For this purpose one can use the decomposition Z^,^ = ip^i, + 61 fj,^. The 
first piece of this decomposition - t/^^,^ is given by Eq. (jH^ . and the second piece - 61 is 
equal to the retarded solution of Eq. (jHTjl . The piece 51^,^ is further decomposed according 
to 61^1, = 6F^^ + 61^1}. Eqs. ()59I61I62I63I64|) provides a prescription for the calculation of the 
desired sufficiently regular piece 61^^ that should be substituted in Eq. (jSHl)- 
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APPENDIX A: LOCAL EXPANSION OF /if J 

In this Appendix we derive Eq. (|13|) . which provides a local expansion of a sufficiently- 
regular tensor- field on a spacelike hypersurface S(r). Most of our derivation follows a similar 
derivation in Ref. We begin our analysis by studying some of the properties of the 
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a{x\x') = \{\i- Ao) / ' g^Mry'^dX , (Al) 



quantities ee*^ and that appear in Eq. (fTSj) . 

First we derive an expression for ee'''^. Consider two points x' and x, and suppose that x' 
hes in a local neighborhood of x, such that within this local neighborhood there is a unique 
geodesic y{\) that connects x' with x. Here A denotes an affine parameter ranging from Ao 
to Al, where x' = y{\o) and x = y{\i)- We denote the square of the invariant length of y{\) 
with 2a{x\x'), and express cr(x|x') as 

^ ;Ai - Anl / 

where y^ = We now calculate ct_^/ = For this purpose we consider an infinitesimal 
displacement of the point x' to a point x' + 5x' , and let ?/(A) + 5y{\) denote the unique 
geodesic connecting x' + 5x' with x. On this geodesic we scale the affine parameter A 
to run from Ao to Ai, such that x' + 5x' = y{Xo) + Sy{\o) and Sy{\i) = 0. Expanding 
6a = a{x\x' + 6x') — a{x\x') to the first order in the variation, and using integration by parts 
we obtain 

6a = (Al - Ao) [g,u{y)yV]',l - (Ai - Ao) ^dX6y>^ + 0{6y') . 

Recalling that y{X) is a geodesic and noting that 6y(Xi) = we obtain 

cr,^,' = -(Al - \o)gi,'u'y''' ■ (A2) 

Hereafter we shall specialize to the space-like geodesic that connects the points z and x on 
S(f). Here we have 2a{y\z) = e'^{y\z). Employing Eq. ()A2|) and using e{y) as the affine 
parameter of the geodesic, we find that 

eVf'e = -e^ , (A3) 
de 

which implies that V'^e is a unit vector, tangent to the geodesic y{e). 

Next we study the properties of the parallel propagator g^. From its definition g^ 
transports an arbitrary vector As, to a vector 

A,iy) = -gJ^iy\z)Asiz), (A4) 

by a parallel propagation of the vector As, on the geodesic y{e), implying that 

\im-gj^[y{e)\z] = 6f,^ . (A5) 



£^0 
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Since -^^Af, 



we find from Eq. ()A4|) tliat 







(A6) 



and g^"g"^ 



The geometric content of implies that g^d^a 

We now construct a local-expansion for h^j^{x) on S(f), this construction follows directly 
from the expansion of h^pijj) on the geodesies y{e). First we introduce the following bi-tensor 
field on y{e) 

Notice that B^^ transforms as a scalar under a coordinate transformation at y. Expanding 
^api^) ~ B&pi^^yi^)] a Taylor series on the geodesic y{e), with z fixed, gives 



de 



+ 0(e' 



J e=0 



Eqs. (IA5IA7jl imply that 



Using Eas. ()A6IA7D we find that 



(A8) 



(A9) 



dB 



a/3 



B,, 



dy^ 



dy^ 



Employing Eqs. ()A3IA5|) we obtain 

-dB 



&f3 



de 



£=0 



-h^:fJz)e'\ 



(AlO) 



We can now rewrite B^^e{%j)\ foTy = x by substituting Eqs. ()A9IA10|) into Eq. ()A8|) . Using 
this expression for B^^ together with Eq. ()A7|) we finally obtain Eq. (fT^. reading 

h'Jix) = -gj^ixlzygfixlz) [h^ - hf^.^^ee^^ + 0{e')_ , 

This expansion takes a very simple form in Fermi-normal coordinates based on zq (For the 
definition and properties of these coordinates see e.g. j^.). In these coordinates we have 
gj' = Sj^ + 0{e^), and Eq. (O is reduced to 

h'J(^) = hlfit) + hllJx^ = 0, t)x^ + 0{e') . (All) 

Here = denotes equality in Fermi normal coordinates, t denotes the time coordinate, x"" 
denote the spatial coordinates, where Latin indices run from 1 to 3. Eq. ()A11|) clearly 
shows that h^^ is a regular tensor-field in a local-neighborhood of the world-line. In Fermi 
normal coordinates it also easy to evaluate the first-order covariant derivatives of Eq. ()A11|) 
since [Vj,= 
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APPENDIX B: LOCAL EXPANSION OF hi 



In this Appendix we study the general form of a local expansion of h^^^ (For more details 
on the construction method of such local expansions see e.g. Ref. 0.)- 

First, we substitute 

Eq. (123) into Eq. (j2HI) and obtain 

V^4x\zGiT)]dT . (Bl) 

Here we introduced the following notation: ar = a{zG{T)\x), ar± = (^)^±, Uf^u[x\zG{T)] = 
U^,uaf3[x\zG{T)]u°'{T)u'^{T), V^^[x\zg{t)] = V^i.a/j [x 1 2:^ (r)]^" (r)u'^ (t) ; the retarded and ad- 
vanced times are denoted r~ and t~^, respectively. These times satisfy cr[x|zG(T''^)] = 0, 
where > f and < f. Introducing the notation Ar^ = — ^|) we expand the 
quantities in Eq. ()B1|) in Taylor series around f and obtain 

a^^ =&^^ a^Ar^ + (l/2)aV(Ar^)2 + 0[{AT^f] , (B2) 
U,,[x\zG{r^)] = U,,[x\z] T Ar^U^,[x\z] + {l/2){Ar^f%,[x\z] + 0[{Ar^f] , 
V^Ax\zg{t)] = V,MA + (r - r)V^A^\z] + (l/2)(r - f)%,[x\z] + 0[{r - rf] . 

From the relation (T[x|2;Gr(r^)] = we find that 

At^ = 6(1 - i^^^,^-^^;^^^^ V) + 0{e') . (B3) 



Notice that expansion ()B3j) was derived using only geometrical consideration. Therefore, 
the higher order terms in this expansion must also have a geometric content. Substituting 
Eq. ()B2|) together with Eq. ()B3|) into Eq. ()B1|) one obtains an expansion in powers of e 
involving only quantities defined on S('r). Each of these quantities can be further expanded 
in powers of e, for example 

In general the coefficients in these expansions are obtained from coincidence limits of the 
expanded quantity and its covariant derivatives . In this example the tensors m^^, m^,^.^, 
and u^^.j^Xf^-j are equal to the coincidence limits (as x z) of Ua/s, Uai3--y, and ?7a/3;(7<5), 
respectively. These coincidence limits have a geometrical content and they are obtained 
from the definitions of a, U^i^a/s and V^i^ais, see e.g. The form of the above mention 
expansions reveals that the expansion of h^^y{x) has the following form 

hl{^)=9r9p'^n=-j'':le-. (B4) 
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Here the expansion coefficients fj^^ are composed only from the following geometrical quan- 
tities: background tensors, m", together with numerical coefficients. The ffist two terms 
in expansion ()B4|) are provided by Eq. (jHUj) . 



APPENDIX C: EXTENDED FERMI GAUGE 

In this Appendix, we provide a prescription for the construction of the extended Fermi 
gauge. In addition, we derive Eq. (jlUj) which provides an expression for h^l^^p in this gauge. 
Our derivation is based on anzats (jH)), reading 



R 

ufia 



;/3 



First, we derive an expression for D^^ap on zcir). Contracting the anzats with and 
comparing with Eq. ()2()|1 we obtain 

D,^^puf' = Q. (CI) 
Taking the covariant derivative V/j of identity (fT^ and using anzats ()44|) we obtain 

Dy[^a]l3 = . (C2) 

Using the identity 

together with anzats (jl^ we obtain 

= ^5<i3 , (C4) 



To define SR^j^^^ we first introduce ^i^R^p-ys which denotes the linear term in the following 
expansion of the Riemann tensor evaluated with the metric g^^ + /i/i^^, (in the Lorenz gauge) 

From which we find that 

Now ^R^l^p is defined by transforming ^R^jj^p Fermi gauge, which gives 

I^^Raf^l = f^i^Ra/B-yS + -^^RafB-ys) = (C5) 

i9>.a + l^h;^i''^)R^f,^s[9,^ + f^h^^i^^] - g.aR'^pJg^.u] + 0(/i2) , 
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R( 

where £^ denotes a Lie derivative, with respect to C,'^. Combining the fact that h^^ = 0, 
together with the fact that h^a^^ satisfies a homogeneous perturbation equation in vacuum, 
reading SR^^ = 0, where 

we find that SR^i^] is traceless. 

Using Eqs. ()C1|C2IC4|) we calculate some of the components of D^^ap- For simplicity we 
express these components in Fermi-normal coordinates, and obtain 

/^.M«t = 0, (C6) 

Here = denotes equality in Fermi normal coordinates. The remaining components are not 
uniquely determined from Eqs. ()(7l|(^2|(74|l . We make the following choice which satisfies 
these equations 

Here spatial components are denoted with Latin indices that run from 1 to 3. To construct 
the extended Fermi gauge one first follows the prescription laid out in Sec. IIB. Next, one 
uses anzats (jl^ together with Eqs. ()C6IC7|) to construct C,u;fiai3- 

Using Eqs. ()(^6p(^7|l together with Eq. (j45|) we obtain Eq. (jl6|) reading 

where 

Here Pj = 5f + u^u^ . 



APPENDIX D: CONSTRUCTION OF 



0/3 



In this Appendix we discuss the solutions to Eqs. (j60|l and (j65|l . First, we consider Eq. 
dnni), reading 
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We substitute Eq. into Eq. (jUU]) . and expand <^^^ in powers of e7l~^, where TZ is the 
smallest of the length scales that characterize Riemann curvature tensor of the background 
geometry. Thus we obtain 

We use Fermi normal coordinates based on the worldline, and expand Qa/s and g"^ in powers 
of e7l~^ as follows 

Notice that the source of Eq. changes on a time scale of 0{TZ). In Fermi normal 

coordinates we have □ = {6°'^dadb) + 0(7^"^). We therefore find that the leading term <^^°j 
satisfies the following Poisson equation 



= e-'A,^ + s-'B,^ + 0,(6°) . (Dl) 
We decompose into spherical harmonics 

oo I 
1=0 m=-l 

and use the spherical-harmonics decompositions We construct the following solution 

for the spherical-harmonics coefficients 

= -[/(/ + l)]-iag + [2 - /(/ + , / > 1 . 

Calculation of a|^^ (using Mathematica software) reveals that it vanishes, and therefore it 



is absent from Eq. ()D3|) . As both sides of Eq. ()D3|) are composed from bi-tensors we 
may replace the notation "=" with "=". Eq. ()D3|) states that 0^ is bounded in a local 
neighborhood of the worldline. Therefore, g^gf^^!^^^ satisfies the original equation ()60|) up 
to an arbitrary 0{e^) term. We therefore set i^^^ = 
We now turn to Eq. (jHSI), which reads 
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Here again we expand {Slap — '^ap) and the corresponding source term in powers of e7l~^. 
Similar to Eq. ()D1|) . the dominant terms as e ^ are obtained by solving the following 
Poisson equation 

(5"U5,)(%-^,^)(°)=p,^. (D4) 

Here dc^gfiSl^^ — <^^^)^°'' is the leading order term in the expansion of {Sla/3 — '^a/s), and 
da'dfPais corresponding leading term in the expansion of the source term, this term 

scales like 5°. We decompose p^^(6', t) and (51^^ — f^a/^Y^^ i^^o spherical harmonics, and 
denote the corresponding spherical-harmonics coefficients with Im indices. Solving for the 
spherical-harmonics coefficients we obtain 

i^h/^ - Vs.p)fl = lil + 3)(2 - ir'piy , / ^ 2 (D5) 
i^ki^ - Va^Y'l = -{l/5)pf^e\l/5 -logs), 1 = 2. 

This solution satisfies the criteria of a sufficiently regular piece. The most general solution to 
Eq. ()D4|) is obtained by adding to Eq. ()D5|) a general homogeneous solution. A homogeneous 
solution which is regular at the worldline has the form of Xll^o Ylm=-i c'™'(^)^/m(6', f)^'', and 
therefore satisfies the criteria of a sufficiently regular piece. Notice that solutions which 
are homogeneous off the worldline but diverge at the worldline are excluded, since they 
correspond to a source in the form of a distribution, and are therefore non- homogeneous. 
The higher-order corrections to {Sla/3 — V&p)^'^'' ^^^'^ sufficiently regular. We therefore 
conclude that {SJa/s — ^ap) satisfies the criteria of a sufficiently regular piece. 
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